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Abstract In this paper we review the possible mechanisms for production of non- 
thermal electrons which are responsible for the observed non-thermal radiation in clus- 
ters of galaxies. Our primary focus is on non-thermal Bremsstrahlung and inverse 
Compton scattering, that produce hard X-ray emission. We first give a brief review 
of acceleration mechanisms and point out that in most astrophysical situations, and 
in particular for the intracluster medium, shocks, turbulence and plasma waves play 
a crucial role. We also outline how the effects of the turbulence can be accounted for. 
Using a generic model for turbulence and acceleration, we then consider two scenarios 
for production of non-thermal radiation. The first is motivated by the possibility that 
hard X-ray emission is due to non-thermal Bremsstrahlung by nonrelativistic particles 
and attempts to produce non-thermal tails by accelerating the electrons from the back- 
ground plasma with an initial Maxwellian distribution. For acceleration rates smaller 
than the Coulomb energy loss rate, the effect of energising the plasma is to primarily 
heat the plasma with little sign of a distinct non-thermal tail. Such tails are discernible 
only for acceleration rates comparable or larger than the Coulomb loss rate. However, 
these tails are accompanied by significant heating and they are present for a short time 
of < 10^ yr, which is also the time that the tail will be thermalised. A longer period of 
acceleration at such rates will result in a runaway situation with most particles being 
accelerated to very high energies. These more exact treatments confirm the difficulty 
with this model, first pointed out by Petrosian (2001). Such non-thermal tails, even if 
possible, can only explain the hard X-ray but not the radio emission which needs GeV 
or higher energy electrons. For these and for production of hard X-rays by the inverse 
Compton model, we need the second scenario where there is injection and subsequent 
acceleration of relativistic electrons. It is shown that a steady state situation, for ex- 
ample arising from secondary electrons produced from cosmic ray proton scattering by 
background protons, will most likely lead to fiatter than required electron spectra or it 
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requires a short escape time of the electrons from the cluster. An episodic injection of 
relativistic electrons, presumably from galaxies or AGN, and/or episodic generation of 
turbulence and shocks by mergers can result in an electron spectrum consistent with 
observations but for only a short period of less than one billion years. 

Keywords intergalactic medium ■ particle acceleration ■ galaxies: clusters: general 



1 Introduction 

In this paper we attempt t o constrain the acce le ration models based o n th e observations 
described in the papers by lOurret et al.l (120081 ) . [llephaeli et al.l (120081 ) and lFerrari et al] 
( 20081 ) - Chapters 4 — 6, t his volume, and the r equired spectrum of the accelerated 
electrons shown in Fig. 6 of IPetrosian et a"l]l2008l - Chapter 10, this volume. However, 
before addressing these details we first compare various acceleration processes and 
stress the importance of plasma waves or turbulence (PWT) as an agent of scattering 
and acceleration, and then describe the basic scenario and equations for treatment 
of these processes. As pointed out below there is growing evidence that PWT plays 
an important role in acceleration of particles in general, and in clusters of galaxies in 
particular. The two most commonly used acceleration mechanisms are the following. 



1.1 Electric field acceleration 

Electric fields parallel to magnetic fields can accelerate charged particles and can arise 
as a result of magnetic field reconnection in a current sheet or other situations. For 
fields less than the so-called Dreicer field, defined as Ej^ = kT/{e\QQ^i), where 

is the collision (electron-electron or proton-proton) mean free patlQ, the rate of ac- 
celeration is less than the rate of collision losses and only a small fraction of the 
particles can be accelerated into a non-thermal tail of energy E < LeEj^. For the ICM 
cm and L ~ 10 cm so that sub-Dreicer fields can only accelerate 
particles up to lOO's of keV, which is far below the lO's of GeV electrons required by 
observations. Super-Dreicer fields, which seem to be present in many sim ulations of 
reconnection (|Drakj . bood : ICassak et all , bood : IZenitani fc Hoshinol boOSi ) , accelerate 



particles at a rate that is faster than the collision or thermalisation time rtijcrm- This 
can lead to a runaway and an unstable electron distribution which, as shown theo- 
retically, by laboratory experim ents and by the above mentione d simulations, most 



retically, by laboratory experim ents and by tne above mentione d 
probably wiU give rise to PWT (|Boris et all Il97d : lHolmanlll985l ). 



In summary the electric Belds arising as a result of reconnection cannot be the sole 
agent of acceleration in the ICM, because there are no large scale magnetically domi- 
nated cosmological Hows, but it may locally produce an unstable particle momentum 
distribution which will produce PWT that can then accelerate particles. 



^ The proton-proton or ion-ion mean free path will be slightly smaller because of the larger 
value of the Coulomb logarithm InTl ~ 40 in the ICM. 
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1.2 Fermi acceleration 

Nowadays this process has been divided into two kinds. In the original Fermi process 
particles of velocity v moving along magnetic field lines (strength B) with a pitch 
angle cos/i undergo random scattering by moving agents with a velocity u. Because 
the head (energy gaining) collisions are more probable than traihng (energy losing) 
collisions, on average, the particles gain energy at a rate proportional to [u/v)^ D^p,, 
where -D/^/j is the pitch angle diffusion rate. This, known as a second order Fermi 
process is what we shaU call stochastic acceleration. In general, the most likely agent 
for scattering is PWT. An alternative process is what is commonly referred to as a 
first order Fermi process, where the actual acceleration occurs when particles cross a 
shock or any region of converging flow. Upon crossing the shock the fractional gain of 
momentum Sp/p oc Ugij/n. Ever since the 1970's, when several authors demonstrated 
that a very simple version of this process leads to a power law spectrum that agrees 
approximately with observations of the cosmic rays, shock acceleration is commonly 
invoked in space and astrophysical plasmas. However, this simple model, though very 
elegant, has some shortcomings specially when applied to electron acceleration in non- 
thermal radiating sources. Moreover, some of the features that make this scenario for 
acceleration of cosmic rays attractive are not present in most radiating sources where 
one needs efficient acceleration of electrons to relativistic energies from a low energy 
reservoir. 

The original, test particle theory of diffusive shock acceler ation (PSA), although 
very elegant and independent of geometry and other details (e.g. iBlandford fc QstrikeJ 
Il978l ) required several conditions such as injection of seed particles and of course tur- 
bulence. A great deal of work has gone into addressing these aspects of the problem 
and there has been a great deal of progress. It is clear that nonlinear effects (see e.g . 
lDrurvlll983l : IBlandford fc Eichle J [l987l : I Jones fc Ellisonlll99ll : iMalkov fc DrurvlboOll ) 



and losses (specially for electrons) play an important role and modify the resultant 
spectra and efficiency of acceleration. Another important point is the source of the 
turbulence or the scattering ag ents. A common practice is to assume Bohm diffu- 
sion (see e.g. lEllison et al.l[2005l ). Second order acceleration effects could modify the 
particle spectra accelerated by shocks (see e.g. ISchUckeiser et all Il993l : l&vk ov et al.l 
I2OOO). Although there are indications that turbulence may be generated by the shocks 
and the accelerated particle upstream, many details (e.g. the nature and spectrum of 
the turbulence) need to be addressed more quantitatively. There has been progress 
on the understanding of generation of the magnetic field and turbulenc e on strong 
shocks (jBell fc Luceal200ll : lAmato fc Blasil . bood : Fviadimirov et"aLl.l2q06l) as required 
in recent observations of supernova remnants (see e.g. Volk et al. 2005h . There is also 



some evidence for these processes fro m observations of heliospheric shocks (see e.g. 
iKennel et al] 1 19861 : lEllison et al.lll99(tl. Basic fe atures of particle acceleration by cos- 



mological shocks were discussed bv lBvkov et al. 1 12008a - Chapter 7, this volume, so we 



will concentrate here on the stochastic acceleration perspective. 



1.3 Stochastic acceleration 

The PWT needed for scattering can also accelerate particles stochastically with a rate 
Dee/E'^i where Dee is the energy diffusion coefficient, so that shocks may not be 
always necessary. In low beta plasmas, /3p = 2{vs/vp^)'^ < 1, where the Alfven velocity 
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v\ — ^ B'^ /A-Kp, the sound velocity Vs — yjYT jm, p = nra is the mass density and 
n is the number density of the gas, and for relativistic particles the PWT-particle 
interactions are dominated by Alfvenic turbulence, in which case the rate of energy 
gain Dee/E^ = {'"a/'")'^ <g ^mmi ^'^ that the first order Fermi process is more 
efficient. However, at low energies and/or in very strongly magnetised plasmas, where 
ua can exceed c, the speed of lighln , the acceleration rate may exceed the scattering 
rate fsee lPrvadko fc Pctrosian'f 997*1. in which case low energy electrons are accelerated 
more efficiently by PWT than by shocks0 

Irrespective of which process dominates the particle acceleration, it is clear that 
PWT has a role in all of them. Thus, understanding of the production of PWT and its 
interaction with particles is extremely important. Moreover, turbulence is expected to 
be present in most astrophysical plasmas including the ICM and in and around merger 
or accretion shocks, because the ordinary and magnetic Reynolds numbers are large. 
Indeed turbulence may be the most efficient channel of energy dissipation. In recent 
years there has been a substantial progress in the understand i ng of MHD turbulence 
dGoldreich fc Sridhaij.[l995l . ll997l : lLithwick fc Goldreichl[20o3 : Icho fc Lazarianl . bood . 
bood l. These provide new tools for a more quantitative investigation of turbulence and 
the role it plays in many astrophysical sources. 

2 Turbulence and stochastic acceleration 

2.1 Basic scenario 

The complete picture of stochastic acceleration by PWT is a complex and not yet fully 
understood or developed process. However, one might envision the following scenario. 

Turbulence or plasma waves can be generated in the ICM on some macroscopic 
scale L ~ 300 kpc (some fraction of the cluster size or some multiple of galactic sizes) 
as a result of merger events or by accretion or merger shocks. That these kind of 
motions or flows with velocity comparable to or somewhat greater than the virial ve- 
locity ~ 1000 kms~^ will lead to PWT is very likely, because in the ICM the 
ordinary Reynolds number Rc = ui^L/v S> 1. Here v ~ ^th-^scat/S is the viscosity, 
«th — -y/kr/m ~ it/,(r/10^)^^^ and Ascat is the mean free path length. The main 
uncertainty here is in the value of Ascat- For Coulomb collisions Agcat ~ 15 kpc (Eq. [TJ 
and Re ~ 1 00 is just barely large enough for generation of turbulence. However, in a 
recent paper lSrunetti fc Lazariar] (2007) argue that in the presence of a magnetic field 
of _B ~ ^G, «A ~ 70(i3//iG)(10~'^ cm~'^/n)^/^ kms~^ is much smaller than ujjj so 
that the turbulence will be super- Alfvenic, in which case the mean free path may be 
two orders of magnitude smallejfl yielding _Rc ~ 10^. We know this also to be true from 
a phenomenological consideration. In a cluster the hot gas is confined by the gravita- 
tional field of the total (dark and 'visible') matter. Relativistic particles, on the other 
hand, can cross the cluster of radius i? on a timescale of Tcross = 3 x 10^ (i?/Mpc) yr 
and can escape the cluster (see Fig. |4] below), unless confined by a chaotic magnetic 
field or a scattering agent such as turbulence with a mean free path Ascat ^ R- If so, 

^ Note that the Alfven group velocity tig = Cy/v^jjv^ + c?) is always less than c. 

^ In pract ice, i.e. mathematically, there is little difference between the two mechanisms 
ijjonesl . fl994l V and the acceleration by turbulence and shocks can be combined (see below). 

* Plasma instabilities, possibly induced by the relativistic particle s , can be another agent of 
decreasing the effective particle mean free path llSchekochihin et al.L l2005f) . 
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then the escape time Tcsc ~ 7cioss(-R/Ascat) = rcross/fscat- The curve marked with 
arrows in this figure shows the maximum value of the required Tgcat so that the es- 
cape time is longer than the energy loss time tiqss- As is evident from this figure, for 
a GeV electron to be confined for a Hubble timescale, or Tesc ~ 10^'' yr, we need 
Tscat ~ 3 X 10*^ yr or Ascat < 10 kpc. This could be the case in a chaotic magnetic field 
an d /or in the presence of turbulence. Some observations r elated to this a. r e disc ussed 
bv lPetrosian et ahlboosl - Chapter 10, this volume; see also IVogt fc Enfilh] l|2005l ). Nu- 



merous numerical simulations also agree with this general picture . There is evidence for 
large scale bulk fiows in the simulations of merging clusters (e.g. IRoettiger et al.lll996l : 
iRicker fc SarazinI l200lh . and that these are converted into turbulence with energies 
that are a substantial fr action of the thermal energy of the clusters ( e.g. [ Sunyacv e t al.l 
I2OO3I : iDolag et al.lboosi '). For more details see lBrunetti fc LazarianI ||2007F . 

Once the PWT is generated it can undergo two kind of interactions. The first is 
dissipationless cascade from wave vectors fcmin ~ to smaller scales. The cascade is 
gouverned by the rates of wave- wave interactions. For example, in the case of weak tur- 
bulence, that can be considered as a superposition of weakly interacting wave packets, 
the three wave interactions can be represented as 

a;(ki) +a;(k2) =a;(k3) and ki+k2=k3, (2) 

where k is the wave vector, and the wave frequency, (jj(k), is obtained from the plasma 
dispersion relation. One can interpret Eq. [2]as energy- momentum conservation laws for 
weakly coupled plasma waves in a close analogy to the optical waves. The interaction 
rates can be represented by the wave diffusion coefficient D^j or the cascade time 
Teas ~ k^/Dij. The largest uncertainty is in the diffusion coefficient. Because of the 
nonlinear nature of the interactions this coefficient depends on the wave spectrum 
W(k). As mentioned above there has been considerable progress in this area in the 
past two decades and there are some recipes how to calculate the diffusion coefficients. 

The second is damping of the PWT by wave-particie interaction which termi- 
nates the dissipationless cascade, say at an outer scale fcmax when the damping rate 
-r(fcmax) = Tclis(fcmax). The range fcmin < k < fcmax is caUed the inertial range. The 
damping rate can be obtained from the finite temperature dispersion relations (see 
below). The energy lost from PWT goes into heating the background plasma and/or 
accelerating particles into a non-thermal tail. These processes are described by the 
diffusion coefficients D^^ and -D/i/i introduced above. These coefficients are obtained 
from consideration of the wave-particle interactions which are often dominated by res- 
onant interactions, specially for low beta (magnetically dominated) plasma, such that 

uj(k) — k cos Ovfj, = nO/^f, (3) 

for waves propagating at an angle 9 with respect to the large scale magnetic field, and a 
particle of velocity v, Lorentz factor 7, pitch angle cos ^ and gyrofrequency Q — eB/mc. 
Both cyclotron (the term in the right hand side of Eq. [3]) and Cerenkov resonance (the 
sec ond term in the left h and side) play important roles in the analysis (see for details 
e.g. lAkhiezer et"al]|l975l ) . Here, when the harmonic number n (not to be confused with 
the density) is equal to zero, the process is referred to as the transit time damping. 
For gyroresonance damping by waves propagating parallel to the field lines {6 = 0) 
n = ±1. For obliquely propagating waves, in principle one gets contributions from all 
harmonics n — ±1, ±2, . . ., but for pra ctical purposes most of the contribution comes 
from the lowest harmonics n = ±1 (see lPrvadko fc Petrosianlll998l ). 
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2.2 Dispersion relations 

It is clear from the above description that at the core of the evaluation of wave-wave 
or wave-particle interactions (and all the coefficients of the kinetic equations described 
below) lies the plasma dispersion relation (jj{k)- It describes the characteristics of the 
waves that can be excited in the plasma, and the rates of wave- wave and wave-particle 
interactions. 

In the MHD regime for a cold plasma 

uj = vp^kcos6 and lj = Vfi^k (4) 

for the Alfven and the fast (magneto-sonic) waves, respectively. Beyond the MHD 
regime a multiplicity of wave modes can be present and the disp ersion relation i s more 
complex and is obtained from the following expressions (see e.g. ISturrock|[l993 ): 



2. ^Pjnj - R){ni - L) 
tan 6 — - — T. 77^5 r, (5) 



where rir = kc/uj \s, the refractive index, S = ^{R~[- L)^ and 

i i i 

(6) 

Here cjpj = 4nniqf/mi and f2i = |gj|_B/mjC are the plasma and gyro frequencies, 
~ <li/\<li\j BtXid Hi, Qi, and rrii are the density, charge, and mass of the background 
particles. For fully ionised plasmas such as that in the ICM it is sufficient to include 
terms due to electron, proton and a particles. Fig. [T] shows the dispersion surfaces 
(depicted by the curves) obtained from the above expressions along with the resonant 
planes in the {lu, fcy, kj_) space. Intersections between the dispersion surfaces and the 
resonant planes define the resonant wave-particle interactions and the particle kinetic 
equation coefficients. One can also envision a similar graphic description of the three 
wave interactions (Eq. [2]) using the intersections of the curved dispersion surfaces. 
However, such calculations have been carried out only in the MHD regime using the 
simp l e relations of Eq. 2] which is already a complicated procedure (see e.g. IChandranI 
I2OO5I : ILuo fc Melrose 2006 ) . 

The above dispersion relations are good approximations for low beta plasmas but 
in the ICM the plasma beta is large: 

/3p = SirnkT/B'^ = 3.4 x 102(n/10"^ cm-^){fiG/Bf{T/W^ K). (7) 

For high beta plasmas the dispersion relation is modified, specially for higher frequen- 
cies UJ ~ fct'th- Foi^ example, in the MHD regime, in addition to the Alfven mode one 
gets fast and slow modes with the dispersion relation (see e.g. ISturrockiil994i ) 



{vl + v^) ± ^vX + vi ^2vlvl cos 261 



(8) 



and the more general di spersion relatio n (Eq. [5j is modified in a more complicated way 
(see e.g. lAndre|[l985l or lSwansonlll989l ). The finite temperature imparts an imaginary 
part L0{ to the wave frequency that gives the damping rate r{k) as long as ui[ < Ur, the 
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Fig. 1 Dispersion relation (curves) surfaces for a cold fully ionised H and He (10 % by num- 
ber) plasma and resonance condition (flat) surfaces showing the regions around the electron 
(top panel) and proton (bottom panel) gyro- frequencies. Only waves with positive 
(or < 6 < tt/2) are shown. The mirror image with respect to the (lu, k±) plane gives 
the waves propagating in the opposite direction. From high to low frequencies, we have one 
of the electromagnetic branches (green), upper- hybrid branch (purple), lower-hybrid branch, 
which also includes the whistler waves (pink), fast- wave branches (yellow), and Alfven branch 
(black). The effects of a finite temperature modify these curves at frequencies u! ~ fc'Uthi where 
titli = ^ 2kT/m is the thermal velocity (see e.g. lAndr3ll985h . The resonance surfaces are for 
electrons with v = 0.3c and |^| = 1.0 (top panel: upper, brown n = 1, lower, light blue n = 0) 
and ■'He (bottom panel: middle, brown n = 1) and ^He (bottom panel: upper, brown n = 1) 
ions with \fj,\ = 1.0 and v = 0.01c. The resonance surfaces for the latter two are the same when 
n = (bottom panel: lower). 
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real p a rt of the frequenc>H. F or more detai ls see e.R. Barnes_&, ScargldJipTal):^ 



1989| ) : IPrvadko fc Petrosiai](|l998l . [l^ Van ballegooiieni (|200a ): lBrunetti fc LazarianI 



20071 ) . In general, these rates and the modification of the dispersion relation are 



known for Maxwellian (sometimes anisotropic) energy distributions of the plasma par- 
ticles. For non-therma l distributions the damping rates can be evaluated as described 
IPetrosian et alJ ( 2006h using the coupling described in Eg. ITT] below. 



2.3 Kinetic equations and their coefficients 
2.3.1 Wave equation 

Adopting the diffusion approximation (see e.g. IZhou fc Matthaeusiri990h . one can ob- 
tain the evolution of the spatially integrated wave spectrum W{]s.,t) from the general 
equation 



dW _ d 
dt dkj 



where Cf^ is the rate of generation of PWT at fcmin, TlSc is the escape time, and 
Dij and F describe the cascade and damping of the waves. The calculation of the 
damping rate is complicated but as described above it is well understood, but there 
are many uncertainties about the treatment of the cascade process or the form of Dij . 
This is primarily because of incompleteness of the theoretical models and sufficient 
observation al or experimental data. There are some direct observations in the Solar 



wind (e .g. iLeamon et al. 199^ ) and indirect inferenc es in the interstellar medium ( 



see e.g. [Armstrong et all I1995I ). There is some hope (jlnogamov fc Sunvaevl l2003h of 



future observations in the ICM. Attempts in fitting th e Solar wind data have provided 
some clues about the cascade diffusion coefficients (see lLeamon et cil]|l999l : |jiang et al.l 
120071 '). 

2.3.2 Particle acceleration and transport 

As described bv lPetrosian et ahlbOOSl - Chapter 10, this volume, the general equation 
for treatment of particles is the Fokker-Planck equation which for ICM conditions can 
be simplified considerably. As pointed out above we expect a short mean free path and 
fast scatterings for all particles. When the scattering time Tgcat = Ascat/t' ~ 0-/D^^j) 
is much less than the dynamic and other timescales, the particles will have an isotropic 
pitch angle distribution. The pitch-angle averaged and spatially integrated particle 
distribution is obtained froir0 



dN{E,t) _ d_ 
dt ^ dE 



DEE ^^f/^ -{A-EL)N{E,t) 



(10) 



Here Dee/E^ is the energy diffusion, due to scattering by PWT a s described above 
and due to Coulomb collisions as discussed bv lPetrosian et ahlbOOSl - Chapter 10, this 



^ Note that the 'thermal' effects change oJr only slightly so that often the real part, the 
resonant interaction rate and the particle diffusion coefficients can be evaluated using the 
simpler cold plasma dispersion relation depicted in Fig. [T] 

^ The derivation of this equation for the stated conditions and some other details can be 
found in the Appendix. 
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volume, A(E)/E ~ i^Dee/E'^, with C(-B) = (2 - 7"^)/(l + 7"^) is the rate of direct 
acceleration due to interactions with PWT and all other agents, e.g., direct first order 
Fermi acceleration by shocks, Ej^jE is the energy loss rate of the particles (due to 
Coulomb collisions and synchrotron and IC losses, see Fig. l4l in IPetrosian et ahlboosl 
- Chapter 10, this volume), and and the term with the escape times TJ^c describe 
the source and leakage of particle^. 

The above two kinetic equations are coupled by the fact that the coefficients of 
one depend on the spectral distribution of the other; the damping rate of the waves 
depends on N{E, t) and the diffusion and accelerations rates of particles depend on 
the wave spectrum W(\i.,t). Conservation of energy requires that the energy lost by 
the waves VVtot = / -r(k)M^(k)d'^fc must be equal to the energy gained by the particles 
from the waves; £ = J[A(E) - As^]N{E)dE. Representing the energy transfer rate 
between the waves and particles by S(k, E) this equality implies that 



r(k)= dEN{E)i:{'k,E), A(E)= d-^kW(k)i:{'k,E) + Ash, (11) 



where we have added Ag^i to represent contributions of other (non-stochastic accelera- 
tion) processes affecting the direct acceleration, e.g., shocks. 

If the damping due to non-thermal particles is important then the wave and particle 
kinetic equations (|9} and (I10|l ar e coupled and attempts have been made to obtain so- 
lutions of the coupled equations ( Miller et al.l . [l99il : lBrunetti fc Blasil . bOOsh . However, 



most often the damping rate is dominated by the background thermal particles so that 
the wave and non-thermal particle kinetic equations decouple. This is a good approx- 
imation in the ICM when dealing with relativistic electrons so that for determination 
of the particle spectra all we need is the boundaries of the inertial range (fcuiin, fcmax), 
the wave spectral index q in this range (most likely 5/3 < q < 3/2), and the shape of 
the spectrum above femax which is somewhat uncertain (see I Jiang et ahlboOTl ). 



3 Particle acceleration in clusters of galaxies 

We now address the problem of particle acceleration in clusters of galaxies. The current 
information on the ICM does not allow us to treat the problem as outlined above by 
solving the coupled kinetic equations. In what follows we make reasonable assumptions 
about the turbulence and the particle diffusion coefficients, and then solve the particle 
kinetic equation to determine N{E, t). We first consider the apparently simple scenario 
of acceleration o f the background thermal particles. Based on some general arguments, 
IPetrosianI 1 200 ll . POl hereafter) showed that this is not a viable mechanism. Here we 



carry out a more accurate calculation and show that this indeed is the case. This leads 
us to consider the transport and acceleration of high energy particles injected into the 
ICM by other processes. 



in what follows we will assume that the waves are confined to the ICM so that Tj^ 
and in some cases we will assume no escape of particles and let T^^c ~> 00 ■ 
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3.1 Acceleration of background particles 

The source particles to be accelerated are the ICM hot electrons subject to diffu- 
sion in energy space by turbulence and Coulomb collisions, acceleration by turbulence 
or shocks, and energy losses due to Coulomb collisionijf]. We start with an ICM of 
kT = 8 keV, n — 10~'^ cm^'^ and assume a continuous injection of turbulence so 
that its density remains constant resulting in a time independent d iffusion and accel- 
eratio n rate. The results described below is from a recent paper bv lPetrosian fc EastI 
( 20071 . '^E07 hereafter). Following this paper we assume a simple but generic energy 
dependence of these coefficients. Specifically we assume a simple acceleration rate or 
timescale 

Tac = E/A{E) = (.D{E)/E^ = ro(l + Ec/Ef. (12) 

Fig.[2]shows a few examples of these time scales along with the effective Coul omb (plus 
IC and synchrotron) loss times as described in Fig. 3 of lPetrosian et al]l2008l - Chapter 
10, this volume. 

We then use Eq. [10] to obtain the time evolution of the particle spectra. After 
each time st ep we use the resulta nt spectrum to update the Coulomb coefficients as 
described bv lPetrosian et aPboOSl - Chapter 10, this volume. At each step the electron 
spectrum can be divided into a quasi-thermal and a 'non-thermal' component. A best 
fit Maxwellian distribution to the quasi-thermal part is obtained, and we determine a 
temperature and the fraction of the thermal electrons. The remainder is labelled as the 
non-thermal tail. (For more details see PE07). The left and middle panels of Fig.[3]show 
two spectral evolutions for two different values of acceleration time to/tcouI — 0.013 
and 2.4, respectively, and for Ec = 25 keV and p = 1. The last spectrum in each case 
is for time t = tq, corresponding to an equal energy input for all cases. The initial 
and final temperatures, the fraction of particles in the quasi-thermal component N^h, 
and the ratio of non-thermal to thermal energies -R„onth are shown for each panel. 
The general feature of these results is that the turbulence causes both acceleration and 
heating in the sense that the spectra at low energies resemble a thermal distribution but 
also have a substantial deviation from this quasi-thermal distribution at high energies 
which can be fitted by a power law over a finite energy range. The distribution is broad 
and continuous, and as time progresses it becomes broader and shifts to higher energies; 
the temperature increases and the non-thermal 'tail' becomes more prominent. There 
is very little of a non-thermal tail for tq > rpoul S'Hd most of the turbulent energy 
goes into heating (middle panel). Note that this also means that for a steady state 
case where the rate of energy gained from turbulence is equal to radiative energy loss 
rate (in this case thermal Bremsstrahlung, with time scale S> tcoui) there will be an 
insignificant non-thermal component. There is no distinct non-thermal tail except at 
unreasonably high acceleration rate (left panel). Even here there is significant heating 
(almost doubling of the temperature) within a short time (~ 3 x 10^ yr). At such rates 
of acceleration most particles will end up at energies much larger than the initial kT 
and in a broad non-thermal distribution. We have also calculated spectra for different 
values of the cutoff energy Ec and index p. As expected for larger (smaller) values of Ec 
and smaller (higher) values of p the fraction of non-thermal particles is lower (higher). 



* In our numerical results we do include synchrotron, IC and Bremsstrahlung losses. But 
these have an insignificant effect in the case of nonrelativistic electrons under investigation 
here. 
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Fig. 2 Acceleration and loss timescales for ICM conditio ns based on the mode l described in 
the text. We use the effective Coulomb loss rate given bv lPetrosian et al 1 l2008l - Chapter 10, 
this volume, and the IC plus synchrotron losses for a CMB temperature of Tcmb = 3 K and 
an ICM magnetic field of B = 1 fiG. We also use the simple acceleration scenario of Eg. ll2l for 
Ec = 0.2meC^ 100 keV) and for the three specified values of p and times ro/rcoul (from 
iPetrosian fc Easdl2007l ). 

The evolution in time of the temperature (in units of its initial value) , the fraction 
of the electrons in the 'non-thermal' component, the energy ratio -Rnonth well as an 
index S = ~d lnAi'(ii')/d IniJ for the non-thermal component are shown in the right 
panel of Fig. [3] All the characteristics described above are more clearly evident in 
this panel and similar ones for p — —1 and +1. In all cases the temperature increases 
by more than a factor of 2. This factor is smaller at higher rates of acceleration. In 
addition, high acceleration rates produce flatter non-thermal tails (smaller 5) and a 
larger fraction of non-thermal particles (smaller N^y^) and energy (-Rnonth)- 

It should be noted that the general aspects of the above behaviour are dictated 
by the Coulomb collisions and are fairly insensitive to the details of the acceleration 
mechanism which can affect the spectral evolution somewhat quantitatively but not 
its qualitative aspects. At low acceleration rates one gets mainly heating and at high 
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Fig. 3 Upper left panel: Evolution witli time of electron spectra in the presence of a constant 
level of turbulence that accelerates electrons according to Eq. 1121 with tq/tq^^i = 0.013, 
Ec = 0.2 (~ 100 keV) and p = 1. For the last spectrum obtained for time t = tq, the low 
end of the spectrum is fitted to a thermal component (dashed curve). The residual 'non- 
thermal' part is shown by the dotted curve. We also give the initial and final values of the 
temperature, the fraction of electrons in the thermal component A^th, and the ratio of energy 
of the non-thermal component to the thermal components i?nonth- Upper right panel: Same 
as above except for tq/tq^^^i = 2.4. Note that now there is only heating and not much of 
acceleration. Lower panel: Evolution with time (in units of tq) of electron spectral parameters, 
T(t)/To, A^tii, ^nonth ^nd the power-law index 5 for indicated values of tq/tcouI a^nd for p = 
and Ec = 100 keV. Note that for models with the sa me value of p at t = tq r oughly the same 
amount of energy has been input into the ICM (from IPetrosian &: Eastil2007r ). 



acceleration rate a prominent non-thermal tail is present but there is also substantial 
heating within one acceleration timescale which for such cases is very short. Clearly in 
a steady state situation there will be an insignificant non-thermal component. These 
findings support qualitatively findings by POl an d do not support the p resence of 
distinct non-thermal tails advocated bv lBlasil ( 2000h and lPogiel et al.l (|2007l ). but agree 
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qualitatively with the more rigorous analysis of IWolfe fc Melial (|2006l ). For further 
results, discussions and comparison with earlier works see PE07. 

We therefore conclude that the acceleration of background electrons stochastically 
or otherwise and non-thermal bremsstrahlung are not a viable mechanism for produc- 
tion of non-thermal hard X-ray excesses observed in some clusters of galaxies. 



3.2 Acceleration of injected particles 

The natural way to overcome the above difficulties is to assume that the radio and 
the hard X-ray radiation are produced by relativistic electrons injected in the ICM, 
the first via synchrotron and the second via the inverse Compton scattering of CMB 
photons. The energy loss rate of relativistic electrons can be approximated by (see POl) 

EL{E)/Ep = {l + {E/Epf)/n,,„ (13) 

where 

Tioss = -Ep/(47rromoc^nln^) and £'p ~ mcC^ [-^ "/"p^ /e ]^^^ (^4) 

axe twice the loss time and the energy where the total loss curve reaches its maximurrj^ 
(see Fig. mi. Here ro = e^/(moC^) = 2.82 x 10~^^ cm is the classical electron radius, 
UpYi (due to the CMB) and B^/Stt are photon (primarily CMB) and magnetic field 
energy densities. For the ICM B ~ fiG, n = 10"'^ cm~'^ and the Coulomb logarithm 
In A = 40 so that tiqss = 6.3 x 10^ yr and Ep = 235mcC^. 

The electrons are scattered and gain energy if there is some turbulence in the ICM. 
The turbulence should be such that it resonates with the injected relativistic electrons 
and not the background thermal nonrelativistic electrons for the reasons described in 
the previous section. Relativistic electrons will interact mainly with low wavevector 
waves in the inertial range where W{k) oc k~'^ with the index q ~ 5/3 or 3/2 for a 
Kolmogorov or Kraichnan cascade. There will be little interaction with nonrelativistic 
background electrons if the turbulence spectrum is cut off above some maximum wave 
vector fcmax whose value depends on viscosity and magnetic field. The coefficients of 
the transport equation (Eq. llOf) can then be approximated by 

D{E) = VE'^,A{E) = aVE''-^, and Tosc = ToscE" . (15) 

For a stochastic acceleration model at relativistic energies a — 2, but if in addition 
to scattering by PWT there are other agents of acceleration (e.g. shocks) then the 
coefficient a will be larger than 2. In this model the escape time is determined by the 
crossing time Tcioss ~ R/c and the scattering time Tgcat ~ ^wl- We can then write 
Tesc ~ 2cross(l + Tcross/Tscat ) . Some examples of these are shown in Fig. |4l However, 
the escape time is also affected by the geometry of the magnetic field (e.g. the degree of 
its entanglement). For this reason we have kept the form of the escape time to be more 
general. In addition to these relations we also need the spectrum and rate of injection 
to obtain the spectrum of radiating electrons. Clearly there are several possibilities. 
We divided it into two categories: steady state and time dependent. In each case we 
first consider only the effects of losses, which means 2? = in the above expressions, 
and then the effects of both acceleration and losses. 



^ We ignore the Bremsstrahlung loss and the weak dependence on E of Coulomb losses at 
nonrelativistic energies. We can also ignore the energy diffusion rate due to Coulomb scattering. 
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Fig. 4 Comparison of the energy dependence of the total loss time (radiative plus Coulomb; 
thick solid line) with timescales for scattering Tscat (lower dashed line), acceleration Tac (thin 
solid lines), crossing time Tcross ^ R/cl3 (long dashed lines) and escape time Tesc ^ Tcrosa /"""scat 
(upper dashed line) of electrons. At low energies the scattering and escape times are dominated 
by Coulomb scattering but at high energies turbulence scattering becomes more important. 
The acceleration is constant at low energies {A{E) oc E) but increases at high energies, cor- 
responding to an acceleration rate A{E) — > a constant at high _E or to a turbulence spectral 
index g = 1. A chaotic magnetic field with scale of 10 kpc will increase the crossing time to 
the primed curve. The arrows show the maximum scattering times for which the escape times 
are equal or longer than the total loss times. The critical energies where the acceleration time 
is equal to the escape time, the Coulomb and inverse Compton loss times are shown. 



3.2.1 Steady state cases 

By steady state we mean variation timescales of order or larger than the Hubble 
time which is also longer than the maximum loss time rioss/2. Given a particle in- 
jection rate Q = Qof(E) (with J f{E)dE — 1) steady state is possible if Tcsc = 
Qo/ J N{E)E-'dE. 

In the absence of acceleration {D — 0) eQ. [10] can be solved analytically. For the 
examples of escape times given in Fig. [4] (Tcsc > tioss) one gets the simple cooling 
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spectra A'' = {Qrioss/ Ep) f{E)dE /{I + (E/Ep)'^), which gives a spectral index 
break at Ep from index po — 1 below to po + 1 above Ep, for an injected power law 
f{E) oc E~P° . For po — 2 this will give a high energy power law in rough agreement 
with the observations but with two caveats. The first is that the spectrum of the 
injected particles must be cutoff below E ~ lOOrricC^ to avoid excessive heating and 
the second is that this scenario cannot produce the broken power law or exponential 
cutoff we need t o explain the radio spectrum of Coma (see Fig. 6 and the discussion in 
IPetrosian et al.ll2008 - Chapter 10, this volume). A break is possible only if the escape 
time is shorter than tq in which case the solution of the kinetic equation for a power 
law injected spectrum (pp > 1 and s > — 1) leads to the broken power law 

Mij?\ n f Tosc(E/Ep)-P°+' HE^Ecr, 
""^""^^^'UossiE/Elr-'-'/is + l) iiE..«E, 



where Ecr = Ep{{s + l)(Tosc/rioss)"^^^^"''^^- Thus, for po ~ 3 and s = and Tc 



CSC 



0.02rjoss we obtain a spec trum with a break at _Ecr ~ 10"*, in agreement with the 



radio data (|Rephaelilll979l model). However, this also means that a large fraction of 



the E < Ep electrons escape from the ICM, or more accurately from the turbulent 
confining region, with a fiux of Fcsc{E) oc N{E)/Tcsc{E). Such a short escape time 
means a scattering time which is only ten times shorter than the crossing time and a 
mean free path of about ~ 0.1-R ~ 100 kpc. This is in disagreement with the Faraday 
rotation observations which imply a tangled magnetic field equivalent to a ten times 
smalle r mean free path. The case for a long escape time was first put forth by I Jaffd 
( 19771 ). 



Thus it appears that in addition to injection of relativistic electrons we also need 
a steady presence or injection of PWT to further scatter and accelerate the electrons. 
The final spectrum of electrons will depend on the acceleration rate and its energy 
dependence. In general, when the acceleration is dominant one expects a power law 
spectrum. Spectral breaks appear at critical energies when this rate becomes equal 
to and smaller than other rates such as the loss or escape rates (see Fig. [Jjl . In the 
energy range where the losses can be ignored electrons injected at energy Eq {f{E) — 
5{E — Eq)) one expects a power law above (and below, which we are not interested 
in) this energy. In the realistic case of long Tosc (and/or when the direct acceleration 
rate is larger than the rate of stochastic acceleration (i.e. a S> 1) then spectral index 
of the electrons will be equal to —q + 1 requiring a turbulen ce spectral index of q = 4 
which is much larger than expected values of 5/3 or 3/2 fsee lPark fc Petrosianlll995l ). 
This spectrum will become steeper (usually cut off exponentially) above the energy 
where the loss time becomes equal to the acceleration time rac = E/A{E) or at Ecr = 
(-EpaDrjoss)^/'-^"'?-' . Steeper spectra below this energy are possible only for shorter 
Tesc- The left panel of Fig. [S] shows the dependence of the spectra on Tcsc for q — 2 
and s = (acceleration and escape times independent of E). The spectral index just 
above Eq is p — y^d/i + 2Tac/2esc — 1.5. In the limit when Tcsc — * oo the distribution 
approaches a relativistic Maxwellian distribution A'^ oc E'^e~^^'^" . For a cut-off energy 
Ecr ~ lO'' this requires an acceleration time of ~ 10^ yr and for a spectral index of 
p = 3 below this energy we need Tcsc ~ Tac/18 ~ 5 x 10^ yr which is comparable to the 
unhindered crossing time. This is too short. As shown in Fig. 2] any scattering mean 
free path (or magnetic field variation scale) less than the cluster size will automatically 
give a longer escape time and a flatter than requir ed spectr um. For fu r ther d etail on 
all aspects of this case see IPark fc PetrosianI l| 19951 ). POl and lLiu et al. I l|2006l ). 
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l°glo('?'/7 ) Electron Kinetic Energy. E (in mr') 

Fig. 5 Left panel: The steady state electron spectra injected at energy Eq and subject to 
continuous acceleration by turbulence with spectral index q = 2 for different values of the 
ratio Tosc/tioss- Note that tq in the label is the same as riass in the text. For very high values 
of this ratio we get a relativistic Maxwellian distribution (from Liu et al. 2006). Right panel: 
Evolution with time of a power law injected spectrum (top line) subject to Coulomb and inverse 
Compton (plus synchrotron) losses as given by Eg. 1171 Solid lines for 6 = 2 (t/n) = 10~"; n = 3, 
2, 1, 0), and dashed Unes for fe = 60 (t/ro = lo-2 i8+"/3. „ = g, 1, 2, 3, 4). The heavy portions 
show the energy range of the electrons needed for production of radio and hard X-rays. 



In summary there are several major difSculties with the steady state model. 



3.2.2 Time dependent models 

We are therefore led to consider time dependent scenarios with time variation shorter 
than the Hubble time. The time dependence may arise from the episodic nature of 
the injection process (e.g. varying A GN activity) and/or from episodic nature of tur- 
bulence generation process (see e.g. ICassano fc Brunettilboosl ). In this case we need 
solutions of the time dependent equation (Eq. I10|> . We start with the generic model 
of a prompt single-epoch injection of electrons with Q{E,t) — Q{E)5{t — to). More 
complex temporal behaviour can be obtained by the convolution of the injection time 
profile with the solutions described below. The results pre sented below ar e from POl. 



Similar treatments o f the f ollowing cases can be found in iBrunetti et all ( 200 ll ) and 
iBrunetti fc LazarianI l|2007l) . 



It is clear that if there is no re-acceleration, electrons will lose energy first at 
highest and lowest energies due to inverse Compton and Coulomb losses, respectively. 
Particles will be peeled away from an initial power law with the low and high energy 
cut-offs moving gradually toward the peak energy Ep. A more varied and complex set 
of spectra can be obtained if we add the effects of diffusion and acceleration. Simple 
analytic solutions for the time dependent case are possible only for special cases. Most 
of the complexity arises because of the diffusion term which plays a v ital role in shaping 
the spectrum for a narrow injection spectrum. For some examples see lPark fc PetrosianI 
Here we limit our discussion to a broad initial electron spectrum in which case 
the effects of this term can be ignored until such features are developed. Thus, if we 
set D{E) = 0, which is a particularly good approximation when a S> 1, and for the 
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purpose of demonstration if we again consider the simple case of constant acceleration 
time {q — 2 and A{E) = aVE), then the solution of Eg. [TOl gives 

cos^ (St / Tiass )[T- (E / Ep) +ta.n(5t/Tiass)/S]P° 

where 5^ = 1 - 6^/4, b = aVroEp = Tiass/rac and r± = 1 ± btan{5t/T^oss)/ i^S)- Note 
that 6 = correspond to the case of no acceleration described above. This solution 
is valid for < 4. For 6^ > 4 we are dealing with an imaginary value for S so that 
tangents and cosines become hyperbolic functions with = 6^/4 — 1. For J = or 
6 = 2 this expression reduces to 

The right panel of Fig. [5] shows the evolution of an initial power law spectrum 
subjected to weak acceleration (6 = 2, solid lines) and a fairly strong rate of acceleration 
(6 = 60, dashed line). As expected with acceleration, one can push the electron spectra 
to higher levels and extend it to higher energies. At low rates of acceleration the 
spectrum evolves toward the generic case of a flat low energy part with a fairly steep 
cutoff above Ep. At higher rates, and for some periods of time comparable to rac, the 
cut off energy Ecr will be greater than Ep and there will be a power law portion below 
it 13 As evident from this figure there are periods of time when in the relevant energy 
range (thick solid lines) the spectra resemble what is needed f or describing the radi o 
and hard X-ray observations from Coma described in Fig. 6 of IPetrosian et ahlbOOSl - 
Chapter 10, this volume. 

In summary, it appears that a steady state model has difficulties and that the most 
likely scenario is episodic injection of relativistic particles and/or turbulence and shocks 
which will re-accelerate the existing or injected relativistic electrons into a spectral 
shape consistent with observations. However these spectra are short lived, lasting for 
periods of less than a billion years. 



4 Summary and conclusion 

We have given a brief overview of particle acceleration in astrophysical plasmas in 
general, and acceleration of electrons in the ICM in particular. We have pointed out 
the crucial role plasma waves and turbulence play in all acceleration mechanisms and 
outlined the equations that describe the generation, cascade and damping of these 
waves and the coupling of these processes to the particle kinetics and energising of the 
plasma and acceleration in both relativistic and nonrelativistic regimes. 

We have applied these ideas to the ICM of clusters of galaxies with the aim of pro- 
duction of electron spectra which can explain the claimed hard X-ray emission either as 
non-thermal bremsstrahlung emission by nonrelativistic electrons or as inverse Comp- 
ton emission via scattering of CMB photons by a population of relativistic electrons. 
It is shown that the first possibility which can come about by accelerating background 
electrons into a non-thermal tail is not a viable mechanism, as was pointed earlier in 

At even later times than shown here on gets a large pile up at the cut off energy (see POl). 
This latter feature is of course ar tificial because we have ne glected the diffusion term which 
will smooth out such features (see iBrunetti fc Lazarianll2007fl . 
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POl. The primary reason for this difhcuhy is due to the short Coulomb coUision and 
loss timescales. Quite generally, it can be stated that at low rates of acceleration one 
obtains a hotter plasma and an insignificant non-thermal tail. Discernible tails can be 
obtained at higher rates of acceleration but only for short periods of time. For periods 
on the order of a billion year such rates will also cause excessive heating and will lead to 
runaway conditions where most of the electrons are accelerated to relativistic energies, 
at which they are no longer bound to the cluster, unless there exists a strong scattering 
agent. 

This leads us to the model where hard X-rays are produced by the inverse Compton 
process and relativistic electrons. Moreover, even if the hard X-ray radiation turns out 
to be not present, or one finds a way to circumvent the above difficulties, we still require 
the presence of relativistic electrons to explain the radio emission. These electrons must 
be injected into the ICM by some other means. They can come from galaxies, specially 
when they are undergoing an active nuclear (or AGN) phase. Or they may be due to 
interactions of cosmic ray protons with thermal protons and the resultant pion decays. 
We have shown that just injection may not be sufficient, because for reasonable injected 
spectra the transport effects in the ICM modify the spectrum such that the effective 
radiating spectrum is inconsistent with what is required. Thus, a re-acceleration in 
the ICM is necessary and turbulence and merger shocks may be the agent of this 
acceleration. In this case, it also appears that a steady state scenario, like the hadronic 
mechanism described above, will in general give a flatter than required spectrum unless 
the electrons escape the ICM unhindered. This requirement is not reasonable because 
the expected tangled magnetic field will increase this time. But, more importantly, the 
presence of turbulence necessary for re-acceleration will result in a short mean free 
path and a much longer escape time. 

A more attractive scenario is if the injection of electrons and/or the production of 
turbulence is episodic. For example for a short lived electron injected phase (from say 
an AGN) but a longer period of presence of turbulence one can determine the spectral 
evolution of the electrons subject to acceleration and losses. We have shown that for 
some periods of time lasting several times the acceleration timescales one can obtain 
electron spectra consistent with what is required by observations. The same will be 
true for a hadronic source if there is a short period of production of turbulence. In 
either case we are dealing with periods on the order of several hundreds of million 
years to a billion years, which is comparable with timescales expected from merging of 
Mpc size clusters with velocities of several thous ands of kms~^ wh ich are theoretically 
reasonable and agree with observations (see e.g. iBradac et al.ll2006l ). 
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and a support from RAS Programs. 



A Particle kinetic equations 

In this section we describe some of the mathematical details required for investigation of the 
acceleration and transport of all charged particles stochastically and by shocks, and the steps 
and conditions that lead to the specific kinetic equations (Eg. llOI I used in this and the previous 
chapters. 
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A.l Stochastic acceleration by turbulence 

In strong magnetic fields, the gyro-radii of particles are much smaller than the scale of the spa- 
tial variation of the field, so that the gyro-phase averaged distribution of the particles depends 
only on four variables: time, spatial coordinate z along the field lines, the momentum p, and 
the pitch angle cos fi. In this case, the evolution of the particle distribution, f{t, z, p, fi), can be 
described by the Fokkor-Planck equation as they undergo stochastic acceleration by interaction 
with plasma turbulence (diffusion coefficients Dpp, D^^ and Dp^), direct acceleration (with 
rate Pg), and suffer losses (with rate pi^) due to other interactions with the plasma particles 
and fields: 



Here /3c is the velocity of the particles and J{t, z,p, fi) is a source term, which could 
be the background plasma or some injected spectrum of particles. The kinetic coefficients 
in the Fokker-Planck equation can be expre s sed through correlation functions of stoc hastic 
electromagnetic fields (see e.g. iMelrosjilQSOl : iBerezinskii et a~ I ll990l : ISchlickeiseHl2002h . The 
effect of the mean magnetic field convergence or divergence can be accounted for by adding 



to the right hand side. 

Pitch-angle isotropy: At high energies and in weakly magnetised plasmas with Alfven 
velocity /3a = va/c <JC 1 the ratio of the energy and pitch angle diffusion rates Dpp/p^ Df^/i ft: 
(/3a //3)^ <SC 1, and one can use the isotropic approximation which leads to the diffusion- 
convection equation (see e.g. iDung &: Petrosianiri994l : iKirk et al]|l988f ): 
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At low energies, as shown by Pryadko & Petrosian (1997), specially for strongly magnetised 
plasmas (a <C 1,/3a > l)i ^vvlv^ 3> f^p, and then stochastic acceleration is more efficient 
than acceleration by shocks (Dpp/ 
be ignored. 



^ . \^,^D^J-f - ^^l.^..n 4- , (23) 

at az op op op 

However. iPetrosian fc Liul ll2004f) find that these dependences are in general weak and one can 
average over the pitch angles. 



A. 2 Acceleration in large scale turbulence and shocks 



In an astrophysical context it often happens that the energy is released at scales much larger 
than the mean free path of energetic particles. If the produced large scale MHD turbulence 
is supersonic and superalfvenic then MHD shocks are present in the system. The particle 
distribution within such a system is highly intermittent. Statistical description of intermittent 
systems differs from the description of homogeneous systems. There are strong fluctuations of 
particle distributio n in shock vicinities . A set of kinetic equations for the intermittent system 
was constructed bv lBvkov fc ToptvginI lll993f) . where the smooth averaged distribution obeys 
an intogro-diffcrential equation (due to strong shocks), and the particle distribution in the 
vicinity of a shock can be calculated once the averaged function was found. 

The pitch-angle averaged distribution function N(r,p,t) of non-thermal particles (with 
energies below some hundreds of GeV range in the cluster case) averaged over an ensemble of 
turbulent motions and shocks satisfies the kinetic equation 



9/ 
dt 



d 



df 



dra or IS 



GLf- 



1 a 



p'^D 



dl 

dp 



+ AL^N + 2BLPf + j(p), 



(24) 



The source term J{t,r,p) is determined by injection of particles. The integro-differential 
operators L and P are given by 



The averaged kinetic coefficients A, B, D, G 
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(25) 



, and Xafi = X ^af) a-re expressed in terms of the 
spectral functions that describe correlations between large scale turbulent motions and shocks, 
the p article spectra index 7 depends on the shock ensemble properties (see lBvkov &: ToptvginI 
Il993f ). The kinetic coefficients satisfy the following rcnormalisation equations; 
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(26) 
(27) 
(28) 
(29) 



Here G = {l/r^h + B). T(k,u!) and S{k,uj) are the transverse and longitudinal parts of the 
Fourier components of the turbulent velocity correlation tensor. Correlations between velocity 
jumps on shock fronts are described by (f){k,ijj), while fi,(k,u)) represents shock-rarefaction cor- 
relations. The introduction of these spectral functions is dictated by the intermittent character 
of a system with shocks. 

The test particle calculations showed that the low energy branch of the particle distri- 
bution would contain a substantial fraction of the free energy of the system after a few ac- 
celeration times. Thus, to calculate the efficiency of the shock turbulence power conversion 
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to the non-thermal particle component, as well as the particle spectra, we have to account 
for th e backreaction of the accelerated particles on the shock turbulence. To do that, iBvko-vj 
1 I2OOID supplied the kinetic equations Egs. I24I I — Il29l with the energy conservation equation for 
the total system including the shock turbulence and the non-thermal particles, resulting in 
temporal evolution of particle spectra. 
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